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ASSESSING MILLISECOND PROTO-MAGNETARS AS 
GRB CENTRAL ENGINES 

Todd A. Thompsoni'2 
RESUMEN 

El resumen sera traducido al espaiiol por los editores. Magnetars are a sizable subclass of the neutron 
star census. Their very high magnetic field strengths are thought to be a consequence of rapid (millisecond) 
rotation at birth in a successful core-collapse supernova. In their first tens of seconds of existence, magnetars 
transition from hot. extended "proto-" magnetars to the cooled and magnetically-dominated objects we identify 
10"' years later as Soft Gamma-ray Repeaters (SGRs) and Anamolous X-ray Pulsars (AXPs). Millisecond 
proto-magnetar winds during this cooling phase likewise transition from non-relativistic and thermally-driven 
to magneto-eentrifugally-driven, and finally to relativistic and Foynting-flux dominated. Here we review the 
basic considerations associated with that transition. In particular, we discuss the spindown of millisecond proto- 
magnetars throughout the Kelvin-Helmholtz cooling epoch. Because of their large reservoir of rotational energy, 
their association with supernovae, and the fact that their winds are expected to become highly relativistic in the 
seconds after their birth, proto-magnetars have been suggested as the central engine of long-duration gamma 
ray bursts. We discuss some of the issues and outstanding questions in assessing them as such. 

ABSTRACT 

Magnetars are a sizable subclass of the neutron star census. Their very high magnetic field strengths are 
thought to be a consequence of rapid (millisecond) rotation at birth in a successful core-collapse supernova. 
In their first tens of seconds of existence, magnetars transition from hot, extended "proto-" magnetars to the 
cooled and magnetically-dominated objects we identify ^ 10^ years later as Soft Gamma-ray Repeaters (SGRs) 
and Anamolous X-ray Pulsars (AXPs). Millisecond proto-magnetar winds during this cooling phase likewise 
transition from non-relativistic and thermally-driven to magneto-eentrifugally-driven, and finally to relativistic 
and Poynting-flux dominated. Here we review the basic considerations associated with that transition. In 
particular, we discuss the spindown of millisecond proto-magnetars throughout the Kelvin-Helmholtz cooling 
epoch. Because of their large reservoir of rotational energy, their association with supernovae, and the fact 
that their winds are expected to become highly relativistic in the seconds after their birth, proto-magnetars 
have been suggested as the central engine of long-duration gamma ray bursts. We discuss some of the issues 
and outstanding questions in assessing them as such. 
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1. INTRODUCTION 

Core-collapse supernovae (SNe) leave behind 
hot proto-neutron stars that cool on the Kelvin- 
Helmholtz timescale (tkh « 10 — 100 s) by radiat- 
ing their gravitational binding energy (~ 10^"^ ergs) 
in neutrinos (Burrows & Lattimer 1986; Pons et 
al. 1999). A fraction of these neutrinos deposit their 
energy in the tenuous and extended atmosphere of 
the PNSH In the standard picture, net neutrino heat- 
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^The charged-current interactions Ueti <-> pe~ and Uep <-> 
ne+ generally dominate heating and cooling. For outflows 
with high entropy per baryon e+e~ <-> up and inelastic 
neutrino-clectron/positron scattering also contribute. 



ing drives a thermal wind that emerges into the post- 
supernova shock environment, blowing a wind-driven 
bubble into the exploding and expanding supernova 
cavity (Woosley et al. 1994; Burrows, Hayes, & Fryx- 
ell 1995). For most massive stellar progenitors with 
extended hydrogen envelopes (Type-II), the cooling 
phase is over well before shock breakout at the pro- 
genitor's surface (~ 1 hour after collapse). For com- 
pact Type-Ibc supernovae, the supernova Shockwave 
traverses the progenitor on a timescale comparable 

to TKH- 

For typical non-rotating non-magnetic (NRNM) 
neutron stars the wind/cooling epoch is in some 
sense a small perturbation to the supernova event 
as a whole. Depending on how one defines the start 
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of the wind phase, the total wind kinetic energy over 
tkh is of order 10*^ — 10^^ ergs, smaU on on the 
scale of the supernova explosion energy, i?gN 10^"'^ 
ergsQ In addition, the total amount of mass ejected 
is a mere ~10~'' — 10""^ M©, a minor addition to 
the few solar masses ejected from the massive stel- 
lar progenitor during the explosion. Finally, because 
of inefficient neutrino heating, the asymptotic wind 
speed does not exceed ~0.1c (Duncan, Shapiro, & 
Wasserman 1986, Qian & Woosley 1996). 

The primary focus of many previous efforts to 
understand proto-neutron star winds — particularly 
in the NRNM limit — has been to assess these out- 
flows as the astrophysical site for production of the 
r-process nuclides (Woosley et al. 1994; Takahashi 
et al. 1994; Qian & Woosley 1996; CardaU & Fuller 
1997; Otsuki et al. 2000; Sumiyoshi et al. 2000; 
Wanajo et al. 2000; Thompson et al. 2001; for a re- 
view, see Thompson 2003a). 

1.1. Rotation & Magnetic Fields 

Simple estimates imply that >10% of all super- 
novae produce magnetars, a class of young neu- 
tron stars (the Anamolous X-ray Pulsars and Soft 
Gamma-ray Repeaters) with inferred large-scale sur- 
face magnetic dipole fields in the range 10^** — 10^^ 
G (Duncan & Thompson 1992; Thompson & Dun- 
can 1993; Kouvehotou et al. 1999; see also Woods & 
Thompson 2004; Kaspi & Helfand 2002). The high 
magnetic fields of magnetars may result from the col- 
lapse of ultra-magnetized iron cores in massive stars 
or white dwarfs in accretion-induced collapse, or they 
may be generated by rapid rotation and efficient dy- 
namo action during the Kelvin-Helmholtz cooling 
epoch of proto-neutron stars (see below. Duncan & 
Thompson 1992). In this latter scenario, it is mil- 
lisecond spin periods at birth that are thought to dis- 
tinguish magnetars from normal neutron stars, with 
their characteristically lower field strengths. 

It is interesting to consider how the standard 
picture of NRNM neutrino-driven protoneutron star 
winds is modified by the presence of a strong large- 
scale magnetic field and rapid rotation. To do so, 
we first quantify the meaning of "strong magnetic 
field" and "rapid rotation" in the proto-neutron star 
context. For this purpose, it useful to consider the 
physical conditions of the proto-neutron star atmo- 
sphere just after birth. 

Assuming a successful explosion occurs, the total 
neutrino luminosity perhaps one second after core 

''However, see Burrows et al. (1995) and Scheck et 
al. (2006) for discussion of how the very early wind is tied 
to the supernova explosion itself. 



collapse is in the range of 5 x 10^^ ergs s~^, di- 
vided roughly equally between neutrino species. The 
temperature at the radius of neutrino decoupling, 
the neutrinosphere, is ^^5 MeV, so that the pressure 
scale height is h fa RNsiRNs/GM){kBT/mp) ~ 0.25 
km for a PNS radius and mass of i?NS — 10 km 
and M ~ 1.4 Mq, respectively. The radius where 
neutrino heating becomes appreciable, the last point 
in the flow where kinetic equilibrium can be main- 
tained, and where the wind-atmosphere becomes ra- 
diation dominated, is characterized by a specific en- 
tropy of s w 4 (see the detailed discussion of Qian & 
Woosley 1996). Assuming that neutrino absorption 
on free nucleons dominates heating, this condition 
determines a characteristic thermal pressure in the 
proto-neutron star atmosphere of 

- 5-3xlO''Lli%,eti%R;,'^' ergs cm-^ (1) 

where we index the total luminosity of the proto- 
neutron star by the luminosity in neutrinos and 
where Lp^,52 = ip,/10^^ ergs s"\ £p,,i5 = ^vj^^ 
MeV is the average energy, and i?io = i?/10 km. 

Although the detailed time evolution of and 
Ey are uncertain, equation ([1]) shows that as the neu- 
tron star cools and Ly and Sy decrease, so too does 
^5=4 Setting B'^ l%Tx equal to yields a crite- 

rion for the magnetic field strength: 

i3,.4 « 9 X 10"L;.f'524f'i5^ro'^' G. (2) 

If i? > -Bs=4, then the magnetic energy density is 
larger than the thermal pressure and it will signif- 
icantly contribute to or dominate the dynamics of 
the outflow. This estimate shows that large-scale 
magnetar-strength magnetic fields of ^^10^^ G are 
required in order to dominate the pressure at very 
early times after explosion. Essential to the following 
discussion, for fixed surface magnetic field strength, 
the magnetic field becomes increasingly dominant as 
the proto-neutron star cools. For example, taking 
Ly oc t^^ and Sy oc L^'*, we find that -85^4 oc t~^\ 
over the few decades of tkh, -83^4 changes apprecia- 
bly. Thus, by "strong magnetic field" we mean in 
this paper that Bns <; ^15=4 (00 

The words "rapid rotation" can be quantified in 
several ways. First, a neutron star with a millisecond 
spin period has a reservoir of rotational energy of 

Srot « 2 X 10^2 ^ p-2 g^gg^ (3) 

^ See, e.g., Pons et al. (1999) for detailed neutrino lumi- 
nosity and average energy time profiles. 

^The s 4 point is generally close enough to the proto- 
neutron star neutrinosphere that it is unnecessary at the level 
of the present discussion to distinguish the surface magnetic 
field strength from Bs=4. 
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where Pi is the spin period in units of 1 ms. There- 
fore, from the standpoint of the evolution of the su- 
pernova remnant, an initial spin period of P < 5 ms 
implies that the rotational energy of the neutron star 
is E^ot ^ EsN ^ 10^^ ergs. That is, if there exists a 
mechanism for spindown of the neutron star on any 
timescale shorter than the remnant age — some way 
to efficiently communicate the energy of rotation to 
the remnant material — a spin period of < 5 ms 
implies that this rotational energy will modify the 
remnant dynamics at order unity or larger. 

A second estimate of what constitutes "rapid 
rotation" , comes from comparing the spin period 
with the characteristic convective eddy turn-over 
timescale (jcon) within the proto-neutron star during 
the cooling epoch. As argued by Duncan & Thomp- 
son (1992), the Rossby number TZ = P/rcon deter- 
mines whether or not an efficient dynamo operates 
in the proto-neutron star; for P ^ 1 ms, the Rossby 
number is less than unity and conditions are ripe for 
dynamo action. Duncan & Thompson (1992) argue 
that it is by virtue of their millisecond rotation at 
birth that magnetars develop 10^"' — 10^^ G magnetic 
fields. 

A last and physically distinct estimate of "rapid 
rotation" comes from noting that the breakup spin 
period of a neutron star is of order 0.5 ms. 

In summary, "strong magnetic field" means B > 
10^^ G at early times after core collapse and the 
criterion for magnetic field domination in the proto- 
neutron star atmosphere decreases as a function of 
time. Additionally, to order of magnitude, "rapid 
rotation" means millisecond spin periods^ An es- 
sential point is that both conditions may be met by 

magnetars at birth. Indeed, Duncan & Thompson 
(1992) argue that P and B are inseparably linked. 

1.2. Millisecond Proto- Magnetars & 
Gamma Ray Bursts 

The combination of rapid rotation and high 
magnetic fields has dramatic consequences for the 
dynamics of any outflow that accompanies proto- 
magnctar cooling. As in studies of magnetic winds 
from rotating stars, the winds from proto-magnetars 
arc dominated by magneto-centrifugal forces. Like 
beads on a wire, the magnetic field lines force the 
wind material into corotation with the stellar sur- 
face out to the Alfven point, where the magnetic en- 
ergy density equals the kinetic energy density of the 

^However, it should be born in mind that both the to- 
tal energy stored in rotation and the properties of the proto- 
neutron star atmosphere as breakup is approached are very 
strong functions of the spin period {P~^ and exponential in 
P~^, respectively). 



outflow. This provides an efficient mechanism for 
spindown (Schatzman 1962; Weber & Davis 1967; 
Mestel 1968ab; Pneuman & Kopp 1971; Belcher & 
Macgregor 1976; Mestel & Spruit 1987). Thus, a 
proto-magnetar's rotational energy can be tapped 
and communicated, in the form of an energetic out- 
flow, to the surrounding medium. 

Although the problem of magnetar birth, and 
how it differs from the birth of more typical neutron 
stars, is interesting in its own right, there are several 
reasons for considering millisecond proto-magnetars 
as the central engine of long-duration GRBs: (1) the 
reservoir of rotation energy is in the range required 
to powers GRBs, (2) proto-magnetar winds, as for 
probably all neutron stars, should become relativis- 
tic on the 10 — 100 second Kelvin-Helmholtz cool- 
ing timescale, and (3) the strong observational con- 
nection between core-collapse supernovae and long- 
duration GRBs is easy to understand in such a 
model. Of course, because magnetars are thought 
to constitute a relatively large fraction of all super- 
novae, not every magnetar can produce a canoni- 
cal long-duration GRB. As in the coUapsar model 
for GRBs (Woosley 1993; Macfadyen & Woosley 
1999), it may be that only millisecond magnetars 
born within compact Type-Ic progenitors produce 
energetic GRBs with ^10^^ ergs. The relativistic 
outflow from millisecond magnetars probably can- 
not sustain high kinetic luminosity for the ~ hours 
required in an extended Typc-II progenitor. Alter- 
natively, it may be that only those magnetars born 
with the highest fields and most rapid rotation yield 
conditions favorable for producing a GRB. 

Assessing millisecond proto-magnetars as GRB 
central engines requires answering a number of ques- 
tions: (1) What is the spindown timescale, and how 
much energy is extracted as a function of time as 
the proto-magnetar cools? (2) How relativistic is 
the outflow and how is the asymptotic Lorentz fac- 
tor connected with the energy loss rate? (3) What is 
the wind geometry, and specifically, how might that 
geometry change as a function of time as the wind 
transitions from non-relativistic to highly Poynting- 
flux dominated? Finally, (4) how are the relativistic 
outflow and the supernova coupled? Are there po- 
tential nucleosynthetic signatures of magnetar spin- 
down? And what might we expect generically from 
the remnants of supernovae that birthed millisec- 
ond proto-magnetars? In this paper we address (1) 
in detail (§2]), drawing on analytic estimates and 
recent exploratory numerical models. We empha- 
size connections with magnetic stellar winds on one 
hand, and force-free models of pulsar spindown on 
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the other — each representing the separate hmits in 
the hfe of a miUisecond magnetar. We discuss (2), 
(3) , and (4) as well as the phases of PNS spindown 
in gl 

2. MILLISECOND PROTO-MAGNETAR 
SPINDOWN 

If the surface magnetic field becomes larger than 
Bs=4, either because L^{t) decreases sufficiently, or 
Bjfs increases sufficiently as a result of dynamo ac- 
tion, the previously solely neutrino-driven outflow 
begins to be accelerated by the action of magneto- 
centrifugal forces. The wind begins to efficiently ex- 
tract rotational energy from the newborn neutron 
star. At this stage the outflow is non-relativistic and 
the timescale for a millisecond magnetar's spin pe- 
riod to e-fold — the spindown timescale — is 



n 
?2 



2 /A/ 
5 I M 



R 



(4) 



where M is the stellar mass, R is the stellar radius. 
M is the mass loss rate, and R^ is the radial position 
of the Alfven point. 

Under the assumption that P « 1 ms, because 
ERot ^ ^^SN (eq- [3]), just one e-folding of ft is suffi- 
cient to modify the dynamics of the supernova rem- 
nant significantly and potentially power a GRB. If rj 
is small with respect to the time for the supernova 
shock to traverse the progenitor (~ tens of seconds 
for type-lb, -Ic progenitors) we also expect this extra 
energy source to modify the supernova nucleosynthe- 
sis (Thompson et al. 2004). 

2.1. The Dominance of Magneto-Centrifugal Forces 

The asymptotic velocity of a wind acceler- 
ated predominantly by magneto-centrifugal forces is 
Voo ~ Ra^- Therefore, a rough criterion for the 
dominance of magneto-centrifugal forces in proto- 
magnetar winds is that Ra^ > V^, where Vi, is the 
asymptotic velocity of an outflow accelerated solely 
by neutrino heating (< 0.1c). For P ^ 1 ms and 
Ra > 15 km this criterion is satisfied, the wind 
is driven primarily by magneto-centrifugal slinging, 
neutrino heating is relatively unimportant in de- 
termining the asymptotic wind velocity, and rota- 
tional energy is transferred efficiently from the proto- 
magnetar to the outflow. 

In order to estimate the spindown timescale and 
the asymptotic velocity of the wind, we must first 
estimate the Alfven radius. In what follows, we dis- 
cuss the non-relativistic limit and then the transition 
to the relativistic regime. We then discuss these es- 
timates in light of recent numerical models, which 
more fully capture the essential physics. 



2.2. The N on- Relativistic Regime: 
Magnetic Stellar Winds 



Angular momentum conservation implies that 
j = d/dt{m) w ~MR\n. The location of the 
Alfven point depends on the radial dependence of 
the poloidal magnetic field. For the purposes of mak- 
ing a simple estimate, we assume that the field is 
strictly monopolar so that B{r) = Bo{RNs/f)^- Us- 
ing p = M/Anr'^Vr and the fact that the Alfven speed 
is VA = Vr{RA) ~ v^{Ra) ~ Ra^, the location of 
Ra is simply 



Ra = B'"^ R^'^Mn)-^'\ 

^ 40 Bli^R%'MZl"Pl'' km. 



(5) 



where Vr is the radial velocity, is the azimuthal 
velocity, p is the mass density, and B is the radial 
magnetic field. In equation ^ Bis = S/10^^ G and 
M is scaled in units of 10^"^ M0 s^^. 

In NRNM protoneutron star winds the mass loss 
rate is (Qian & Woosley 1996) 



Mnrnm « 4 X efX R%^ Mq s" 



(6) 



for a 1.4 Mq neutron star, considerably less than 
that implied by the scaling of equation When 
the magnetic field is strong, centrifugal forces ex- 
tend the scale height of the proto-magnetar atmo- 
sphere. Because the mass outflow rate is essentially 
determined by the density proflle at the sonic point 
[Rs), this effect can increase M significantly if Ra is 
larger than the nominal Rg. For a wind dominated 
by magneto-centrifugal acceleration, one finds that 



Rs w {GM/n^f/^ 



n pI'^ km. 



(7) 



Because equation ([5]) implies that Ra > Rs for fidu- 
cial values of B and M, the larger scaling is in part 
justified. Indeed, more detailed calculations show 
that for P = 1 ms and Lp^,52 = 1, M « 10"^ Mq 
s-i (Thompson et al. 2004; Metzger et al. 2006). 

Using equation ^ we estimate that the absolute 
value of the rotational energy loss rate is 



Pi 



NR. 



54/3^8/3^1/3^4/3 



4/3„8/3,Vl/3p-4/3 e,g3 3-1^(8) 



io5iprr^i^^Af_3 



The subscript 'NR' is added to emphasize that when 
the flow is non-relativistic, P depends explicitly 
on M. The spindown timescale fl/Q in the non- 
relativistic limit is 



r,j^„ ~ 30 s MiaM: 



1/3^-2/3 p-4/3 p-2/3 



10 



P 



15 



p; 



(9) 
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Note that for slower rotation, larger spin periods, the 
spindown tiniescale decreases at fixed M. 

There are several uncertainties in these simple 
estimates. First, the wind is not isotropic, so that 
the factor of Att that appears in the relation between 
M and p is incorrect. Second, the radial velocity, 
even in the idealized problem presented here is not 
Ra^, but can differ at the factor of two level. Most 
importantly, the field is assumed monopolar, when 
in reality the surface field must be dipolar at low- 
est order. However, one may argue that even if the 
star has a dipole field it is the open field lines that 
carry the energy and angular momentum, and per- 
haps these may be — to first approximation — con- 
sidered monopolar. For this reason, it turns out that 
the assumption of a monopole field is remarkably 
good, as long as the surface field strength is suitably 
renormalized to refiect only the field lines that are 
opened to inifinity by the wind. We return to this 
issue below in discussing numerical models. 

2.3. The Non-Relativistic to Relativistic Transition 
As the neutrino luminosity decreases, the char- 
acteristic thermal pressure in the atmosphere of a 
proto-magnetar decreases (eq. [T|). For fixed sur- 
face magnetic field strength, we thus expect mil- 
lisecond magnctar winds to become increasingly 
magnetically-dominated. In particular, M in equa- 
tion dSI should decrease as decreases (see eq. |6|) 
and therefore Ra should increase. Although 
can decrease arbitrarily as < — > oo, Ra cannot in- 
crease indefinitely; instead, it approaches the ra- 
dius of the light cylinder i?L = c/J7 ^ 48Pi km 
asymptotically. As it does so, the fiow becomes 
increasingly relativistic. This is the transition be- 
tween non-relativistic magnetically-dominated mass- 
loaded outflow and relativistic Poynting-flux domi- 
nated wind. All neutron stars likely pass through 
such a transition, regardless of their initial spin pe- 
riod and magnetic field strength. Millisecond proto- 
magnetars are particularly interesting candi- 
date central engine for GRBs because this transition 
to relativistic fiow occurs on the Kelvin-Helmholtz 
timescale (~ 10 — 100 s) and at high wind kinetic 
luminosity, and because magnetars are born in su- 
pernovae. 

Setting Ra = -Rl in equation ^ we can estimate 
the critical mass loss rate below which the wind be- 
comes relativistic: 

Merit = B^R'^n^c-^ 

- 7 X lO-'^Blr,RfoPi^ Mq s"\ (10) 

where B refers to the equivalent surface monopole 
field strength. For lower effective field strengths. 



Merit decreases — that is, the mass flux must 
decrease further in order to enter the relativistic 
regime. For a neutron star born with a 10 ms spin 
period and a 10^^ G surface field. Merit ~ 7 x 10^^^ 
M0 s~^. Thus, for weaker fields, the transition oc- 
curs at lower M, lower L,^, and at a time longer after 
collapse and explosion, a time later in the Kelvin- 
Helmholtz cooling epoch. 

2.4. Relativistic Winds & The Force-Free Limit 

Once Ra becomes close to Rl, the degree of 
Poynting-flux domination is quantified by the pa- 
rameter (Michel 1969; Goldreich & Juhan 1970) 



ct{Rl 



B' 



(11) 



where 7 is the Lorentz factor and B = B{Rl). If the 
outflow is driven primarily by magneto-centrifugal 
forces, ^(Rl) ~ 1. Roughly speaking, if energy 
transfer from the electromagnetic field to the matter 
is efficient, (T(i?i) measures the maximum acheivable 
asymptotic Lorentz factor of the outflow as r — » cxd. 
Assuming that ^{Rl) ^ 1, (^{Rl) can be written 
simply in terms of the mass loss rate: 



(t{Rl) 



B'' R^g^^c 



70 Bf.RtoPr'M: 



(12) 



where M„5 = M/10^^ Mq s~i and we have again 
assumed a monopole field geometry. Compare equa- 
tion with equation pO|) : setting M = Merit in 
the former, yields ij{Rl) = 1. Thus, as L,^ continues 
to decrease as the cooling epoch progresses, we ex- 
pect the flow to become increasingly relativistic and 
Poynting-flux dominated (cr increases). 

In the limiting case of force-free electrodynam- 
ics one neglects the inertia of the matter completely. 
This limit corresponds formally to cr ^ 00. Recent 
numerical models have calculated the spindown lu- 
minosity for a neutron star with aligned magnetic 
and rotational axes in the force-free limit (Con- 
topolous et al. 1999; Gruzinov 2005; Komissarov 
2006; Spitkovsky 2006; McKinney 2006); they flnd 
that 

- 5.8 X W^^Bf^R^P^^ ergs s'^ (13) 

where B here corresponds to the surface dipole field 
strength. In the case of an orthogonal rotator — 
with the magnetic axis perpendicular to the rotation 
axis — Ep-p is larger than that in equation (|13p by 
precisely a factor of two (Spitkovsky 2006). 
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In the force- free limit, the spindown timcscalc for 
the ahgncd rotator is 

rj^^ ~ 760 s Afi.4 i?ro* ^if (14) 

It is tempting to compare E-p-p and rj^p with equa- 
tions ([8|) and Q directly. Unfortunately, such a com- 
parison is complicated by the fact that the former 
are written assuming a dipole surface field strength, 
while the latter assume a monopole field strength. 
One requires knowledge of the amount of open mag- 
netic flux at a given M to make a fair comparison. As 
we discuss in more detail below, the results of Buc- 
ciantini et al. (2006) show that for typical millisecond 
proto-magnetar parameters the equivalent monopole 
field strength is a factor of ~ 3 less than the dipole 
field strength. Therefore, to compare -£/nr and -E'ff ? 
the former should be decreased by a factor of 4, 
implying that E^k/Eyf ~ m^^^'^ m]I^ P^'^ . This 
expression has a strong dependence on P and a 
rather weak dependence on M . For example, tak- 
ing P = 3 ms one finds that the ratio ^^nrZ-Er is 
nearly 20 times larger. Thus, the monopole scal- 
ings derived in ^2.2\ suggest that the early spindown 
of proto-magnetars is rapid because the flow is non- 
relativistic, but magneto-centrifugally driven. The 
spindown luminosity is significantly larger in this 
phase than an application of the force-free spindown 
law would suggest. 

The original suggestion that rapidly rotating, 
highly magnetic neutron stars might power GRBs 
in Usov (1992) essentially employed the force-free 
approximation. Similar, but more fully conceived, 
models were developed by Thompson (1994) and 
Wheeler et al. (2000). The latter addressed the non- 
relativistic wind phase, but with a dipole-likc scal- 
ing for the magnetic field strength. Thompson et 
al. (2004), whose analytic estimates we have so far 
in essence summarized also constructed approximate 
numerical solutions of magneto-centrifugal winds 
from proto-magnetars in order to assess the impor- 
tance of these forces in setting the mass loss rate and 
the spindown luminosity. 

2.5. Numerical Results 

Two sets of results for the spindown luminosity 
have been presented in the preceding discussion: (1) 
the non-relativistic, but magneto-centrifugally domi- 
nated monopole (eq.[8]) and (2) the ahgned force- free 
rotator (eq. [13]). At face value, the latter is more 
secure, coming directly from recent numerical calcu- 
lations. However, the applicability of the force-free 
result is suspect because it is only formally valid in 



the limit a ^ oo. We expect quantitative and poten- 
tially qualitative differences between this limit and 
the case where a is merely, say, 2, 10, or 100. We 
wish to understand the applicability and limitations 
of the two limits in the millisecond proto-magnetar 
context. A set of recent numerical calculations ad- 
dress these issues directly. 

Metzger et al. (2006) solve the time-dependent 
one-dimensional Weber & Davis (1967) problem of a 
non-relativistic magnetically-dominated flow in the 
equatorial plane, including neutrino heating and 
cooling and an appropriate equation of state. This 
approach necessarily employs a monopole magnetic 
field configuration. They provide a detailed account 
of the physics of millisecond proto-magnetar spin- 
down from thermal-magneto-centrifugal winds. We 
emphasize just two results from the study of Met- 
zger et al. (2006): (1) because Ra^ overestimates 
the radial velocity at the Alfven point, the spindown 
timescale in the monopole limit is generally a factor 
of 1.5 — 2 lower than that given in equation Q and 
(2) the mass loss rate increases exponentially when 
Ra is larger than the sonic point (eq. [7]) so that we do 
expect an epoch of enhanced mass loss as the proto- 
magnetar wind becomes increasingly magnetically- 
dominated, even as L,y decreases. 

Bucciantini et al. (2006) (B06) solve the 
two-dimensional axisymmetric time-dependent wind 
problem with both monopole and dipole surface 
magnetic fields, in general relativity. The mass out- 
flow rate is determined self-consistently by imposing 
a flnite thermal pressure at the neutron star sur- 
face, as appropriate in the proto-magnetar context. 
They incude the effects of neutrino heating and cool- 
ing in a parameterized way by employing a F-law 
equation of state. B06 explore both low magnetiza- 
tion non-relativistic winds and relativistic Poynting- 
flux dominated outflows with cr > 1. For the non- 
relativistic, but magnetically-dominated monopole, 
they flnd good agreement with equation ([H]). In the 
high-(T regime, they confirm the analytic force-free 
monopole limit (Michel 1991; Beskin et al. 1998). 

B06 find very interesting results for the aligned 
dipole. For a strong surface dipole field a region 
of closed magnetic field fines (the "closed zone") 
develops at mid-latitudes around the equator in a 
helmet streamer-like configuration. At higher lati- 
tudes, the field lines are opened to infinity by the 
outflow. The last closed fleld line meets the equato- 
rial plane at a y-type point and its radial position 
in that plane is denoted Ry ■ In the force-free limit 
one expects Ry ~ Rl at the equator (e.g., Con- 
topolous et al. 1999). In the dipole case, there is not 
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a onc-to-onc correspondence between the magnitude 
of the surface magnetic field strength and the open 
magnetic flux because as B is increased, the closed 
zone becomes larger. For this reason, B06 find it 
convenient to write an approximate relation, derived 
from the simulations, that relates the dipole surface 
field strength at the pole to an equivalent monopole 
surface field strength: 

Br{RNS,Q = 0) « Br,cq-miRNs) X ^"^^ — ^ ■ 

(15) 

B06 find that the spindown luminosity of the aligned 
dipole with surface field Bj.{Rns,(^ = 0) is the 
same as that for the 2D monopole if they normal- 
ize in terms of the open magnetic fiux. Thus, the 
aligned dipole with Br{RNSiO = 0) spins down hke 
a monopole with Br^cq-miRNs) given by equation 
(fT5|) . Importantly, for a as large as ^ 20, B06 find 
that Ry is considerably less than Rl and that the 
ratio Ry/Rl is a weak function of the magnetic 
field strength. For parameters typical of millisecond 
proto-magnetars, they find that Ry/Rl ^ 1/4 — 1/2. 
The direct implication of this result is that the 
monopole scalings for spindown derived in t j2.2l are 
applicable to spindown with dipole fields, but that 
the magnetic field strength should be decreased in, 
for example, equation ([8]) by a factor given by equa- 
tion (fT5|). Thus, if one assumes a dipole field strength 



of 10 G and that Ry/Rns = 2 (compare with 
B06), then a field strength of should ap- 

pear in equation ([8]). Indeed, BOG find that even as 
(T increases to values as large as ~ 20, spindown is 
more efficient than an application of the pure dipole 
force-free limit would predict (cq. fO)) because Ry 
is significantly less than Rl. These results suggest 
that only at significantly higher a does Ry approach 
Rl and the force-free limit obtain. In the context of 
proto-magnetar spindown, this means that the tran- 
sition to the force- free limit occurs at smaller L^{t), 
later in the Kclvin-Helmholtz cooling epoch. 

3. DISCUSSION & CONCLUSIONS: 
ARE MILLISECOND PROTO-MAGNETARS 
GRB CENTRAL ENGINES? 

3.1. Phases of Spindown Evolution 

There are five phases of spindown in any very 
young rotating neutron star's life (see also B06): (1) 
a pressure-dominated phase in which the wind is 
driven by neutrino-heating {B < Bs=i] eq. H]), (2) a 
phase in which magnetic field effects are present, but 
not dominant so that Ra^ < 0.1c « ct, where ct 
is the isothermal sound speed at the proto-neutron 



star surface, (3) a non-relativistic magnetically- 
dominated phase when Ra is greater than R and 
Rs, but less than R^, (4) a relativistic phase in which 
R A Rl, hut Ry < Rl (as in B06), and lastly (5) 
an epoch when the force-free limit is (presumably) 
applicable and Ry — Ra — Rl- Roughly speak- 
ing, for monotonically decreasing L^, phases (1) — (5) 
represent a time evolution starting immediately after 
the supernova explosion commences. 

The transition from phase (1) to phase (2) is dic- 
tated by equation (HI). When B'^/Stt exceeds Ps=4 
near the proto-magnetar surface, phase (2) begins. 
As the magnetic field becomes increasingly impor- 
tant, the sonic point will move to smaller radii and 
the radial scale over which the magnetic field domi- 
nates will move to larger radii. When Ra is of order 
Rs, phase (3) begins. This transition will be compli- 
cated by the fact that M may increase as a result of 
increased magneto-centrifugal support in the quasi- 
hydrostatic atmosphere (see Metzger et al. 2006 for 
details). Throughout phase (3) Ra^ > ct, < 1, 
and Ra increases from of order the slow magne- 
tosonic radius to Rl- The characteristic spindown 
luminosity is 

E^ix 10^" Bl^^Rt„P~^^^M~l^'' ergs s-\ (16) 

where B is the equivalent monopole field (see ea. [T5)) . 
Because we expect this phase to last of order the 
Kclvin-Helmholtz timescale, the total amount of en- 
ergy extracted is comparable to the asymptotic su- 
pernova energy, ~ 10'"'^ ergs. 

The transition from phase (3) to phase (4) oc- 
curs at (T w 1 {Ra ^ Rl) and for a critical M given 
by equation (fTOjl . Throughout phase (4) the fiow 
is relativistic and Poynting-fiux dominated (cr > 1). 
The total energy and angular momentum loss rates 
are smaller than in phase (3). However, because 
Ry < Rl, the spindown rate is larger than what 
would be inferred from an application of the force- 
free limit. The results of B06 indicate that this is 
true despite the fact that a is larger than and, in- 
creasingly with time as the neutrino luminosity de- 
creases, much larger than unity. Eventually, Li, and 
M decrease sufficiently that Ry ~ Rl, and phase (5) 
begins. A simple and very uncertain extrapolation of 
the results of B06 suggest that this transition occurs 
at a very high a (perhaps ~ 10^). In this epoch, the 
spindown luminosity is given by the recent force-free 
calculations described in tj2.4l 

There is a last phase (or set of phases) beyond 
the scope of the present work, which follow after the 
Kelvin-Helmholtz cooling epoch, as the MHD ap- 
proximation in the magnetospherc of any young ncu- 
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tron star or magnetar begins to break down, the flow 
becomes charge-separated, and particles are acceler- 
ated electromagnetically and to high Lorentz factors 
directly off of the neutron star surface. 

3.2. Energetics, Relativity, & Variability 

As first emphasized by Usov (1992), the energy 
budget and luminosity of millisecond magnetar spin- 
down is in the range needed to explain cosmological 
long-duration GRBs (eqs. [3] & [131 see also Thomp- 
son 1994; Wheeler et al. 2000; similar models also 
by Katz 1997; Kluzniak & Ruderman 1998). The 
recent work of Thompson et al. (2004), Metzger et 
al. (2006), and 606, which we have reviewed here, 
have explored a more complete picture of proto- 
magnetar spindown, from the initial non-relativistic 
wind stage through the Kelvin-Helmholtz cooling 
epoch. 

In order to power a GRB, the outflow must si- 
multaneously achieve high spindown luminosity and 
high Lorentz factor (700 > 100; e.g., Lithwick & 
Sari 2001). The asymptotic Lorentz factor depends 
on the magnetization of the flow, measured by a 
(eq. Ilip . which is, in turn, set by the time depen- 
dence of M (eq. [T^ . itself determined by Li,{t). 
Thus, the transition to a relativistic flow {a > 1) is 
governed by the Kelvin-Helmholtz cooling timescale 
tkh ~ 10 — 100s. If we take the average duration 
of a long GRB to be tgrb ^ 30s, the constraints 
that (1) the total energy ejected must be of order 
£^GRB ^ 10^^-^ ergs and (2) that 7cx, ^ 100 imply 
that there must be efficient conversion of electromag- 
netic energy to kinetic energy in the outflow beyond 
Rl- The reason follows from noting that a oc 
and that E is a. decreasing function of M before the 
force-free limit is reached, after which E is indepen- 
dent of M. These facts limit the space of possible a 
obtainable at a given E to the range of several hun- 
dred, depending on B and at the proto-magnetar 
surface. Therefore, in order to have the Lorentz 
factors required for GRBs, 700 ~ (j{Rl)- This is 
in contradiction to expectations from the force-free 
monopole, which gives 700 ~ g^Rl)^/^ and is poten- 
tially accomplished by efficient magnetic dissipation 
in the out-flowing wind at radii larger than or com- 
parable to Rl and/or the fast magnetosonic point 
(Drcnkhahn & Spruit 2002; Lyutikov & Blandford 
2003). 

We expect variability in both E and a due to 
rapid changes in mass loading. In an average sense, 
this may cause the wind to alternate rapidly between 
(T < 1 and (T > 1. Proto-magnetars, like all neutron 
stars, are expected to be fully convective during tkh 



(Keil et al. 1996; Pons et al. 1999; Thompson & Mur- 
ray 2001; Dessart et al. 2006a). Large-scale convec- 
tion, magnetic field footpoint motion, and instabil- 
ities in the magnetosphere may all cause variability 
in these important quantities [E, 7, a) on millisec- 
ond timescales. Strong variations in the mass load- 
ing could be caused by shearing of large-scale closed 
magnetic loops on the surface of the fully convective 
millisecond proto-magnetar core (Thompson 1994). 
Internal shocks in the flow are a natural consequence 
of these rapid changes in M (Rees & Meszaros 1994). 

3.3. Emergence & Geometry 

Evidence for coUimation abounds in afterglow 
observations of GRBs. On the other hand, there 
is much theoretical work supporting the conclusion 
that it is difficult to coUimate Poynting-flux domi- 
nated flows (e.g., Begelman & Li 1994; Lyubarsky 
& Eichler 2001). If the relativistic proto-magnetar 
wind can be collimated, we expect its emergence 
from the progenitor to resemble models of coUap- 
sar jets escaping their Type-Ibc hosts (e.g., Aloy et 
al. 2000; Zhang et al. 2003, 2004). But, can relativis- 
tic proto-magnetar winds be collimated? 

There are a number of possible answers. The 
first potential answer may be that coUimation sim- 
ply cannot be accounted for in the millisecond proto- 
magnetar model for GRBs. The second reply is that 
relativistic Poynting-flux dominated winds actually 
can be efficiently collimated (see Vlahakis & Konigl 
2003; Vlahakis 2004). A third possibility is that a 
small disk forms outside the rapidly rotating magne- 
tar and that this aids coUimation, or that the proto- 
magnetar is so distorted by centrifugal forces that 
it is disk- like at early times (e.g., see Figs. 7 & 8 
of Dessart et al. 2006b). This option might provide 
a picture which connects logically with the coUap- 
sar model of GRBs (Macfadyen & Woosley 1999). 
A fourth option, as suggested by Thompson (2005), 
is to appeal to a somewhat different geometry. Pul- 
sars drive energetically dominant high Lorentz factor 
equatorial outflows (Komissarov & Lyubarsky 2004; 
Spitkovsky & Arons 2004). The models of B06 show 
that when a > 1, 7(6*) is peaked in the equatorial 
plane. In addition, the presence of the current sheet 
in this region may facilitate the magnetic dissipation 
required for efficient conversion of magnetic energy 
to kinetic energy. In analogy with pulsar winds, per- 
haps it is possible that the geometry of many GRBs 
is "sheet"- or "fan" -like rather than jet-like so that 
the solid angle subtended by the GRB is ^ 9 in- 
stead of ~ 9^. Modeling shows that the expected 
sheet-break (the analog of the jet-break; e.g., Rhoads 
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1999; Sari et al. 1999) is not steep enough to explain 
all achromatic breaks in GRB afterglow lightcurves 
(T. Thompson, unpublished; Granot 2005). In addi- 
tion, for the same observed isotropic equivalent en- 
ergy and break time, a sheet-like geometry increases 
the true energy of the burst with respect to a uni- 
form jet by a factor (x 0'-^. Nevertheless, this geom- 
etry remains an interesting alternative for a subset of 
bursts with shallower break profiles (see, e.g., Fig. 2 
of Panaitescu 2005). 

There is a another interesting possibility that 
may bear on the question of coUimation in mil- 
lisecond proto-magnetar winds. The relativistic 
Poynting-flux dominated wind of phases (4) and (5) 
comes after the non-relativistic mass-loaded a < 1 
wind of phase (3) (see ij3.ip . This means that the 
envelope structure that the relativistic wind encoun- 
ters has been "pre-processed" by the preceding non- 
relativistic flow. The results of BOG (in analogy with 
models of non-relativistic stellar winds; e.g.. Smith 
1998) show that the energetic flux is strongly di- 
rected along the rotation axis by hoop stress when 
a < 1. Thus, if the total energy extracted from 
the proto-magnetar in phase (3) is on the order 
of the supernova explosion energy, then there will 
be a relatively "hollow," asymmetric, and elongated 
channel that the subsequent relativistic flow emerges 
into. Pressure forces and wind material bounding 
this channel may force the less-energetic and rela- 
tivistic flow into a jet-like structure. 

3.4. Remnants & Nucleosynthesis 
For fiducial proto-magnetar parameters, the 
timescale for extraction of an amount of rotational 
energy comparable to 10^^ ergs is small on the 
timescale for the supernova Shockwave to traverse 
a Type-II supernova progenitor and comparable to 
that for a Type-Ibc progenitor. Much of this en- 
ergy is extracted during phase (3) {a < 1), when the 
energetic fiux is strongly directed along the axis of 
rotation because of hoop stress (as in B06) . The ac- 
tion of this outfiow may cause an asymmetry in the 
supernova remnant (see also Wheeler et al. 2000). 
B06 show that the zenith angle at which the ener- 
getic fiux is maximized is an increasing function of 
a so that for a > 1, the energetic loss is primarily 
equatorial. Depending on the timing of the start of 
phase (3) with respect to the position and energy 
of the preceding supernova Shockwave, the action of 
the energetic wind could modify the nucleosynthesis 
in the remnant in an asymmetric way. 

As a bit of speculation and for illustrative pur- 
poses, we note that the Cass A supernova rem- 
nant has a strong jet/counter-jet morphology with 



a distinctive nucleosynthetic enrichment signature 
(Hwang et al. 2001, 2004; WiUingale et al. 2002; Fe- 
scn et al. 2006) and that it has been at least dis- 
cussed in the context of GRB remnants (Laming 
et al. 2006). Indeed, many core-collapse supernova 
events exhibit asymmetry in their spectropolarime- 
try (e.g., Wang et al. 2001). However, for large-scale 
asymmetries representing total energy comparable to 
the asymptotic supernova remnant energetics (as in 
Cass A), we require surface magnetic fields of > lO^"' 
G and rotation rates at birth in the millisecond 
range. Interestingly, Chakrabarty et al. (2001) sug- 
gest that the X-ray point source in the Cass A rem- 
nant is an Anomalous X-ray Pulsar, a magnetar, and 
so perhaps this object satisfied our requirements at 
birth. Evidence also exists for infrared light echoes 
from the X-ray point source, which are interpreted 
as burst-like events, potentially analogous to the gi- 
ant fiares seen from magnetars (Krause et al. 2005; 
as in, e.g., SGR 1806-20, Palmer et al. 2005). 

Because > 10'^^ ergs can be extracted on 
a timescale shorter than or comparable to the 
timescale for the supernova Shockwave to traverse 
the progenitor, we expect that the wind may signif- 
icantly affect the nucleosynthetic yield and its an- 
gular distribution. If significant energy can be ex- 
tracted and communicated to the surrounding en- 
velope of expanding supernova shocked gas rapidly, 
the ^^Ni yield of proto-magnetars may be enhanced. 
In this way it may be possible to generate hyper- 
energetic or 1998bw-like supernovae (Thompson et 
al. 2004; Woosley & Heger 2003). The inferred en- 
ergetics and s^Ni yield of SN2003dh and SN1998bw 
put strong constraints on any GRB mechanism. In 
the coUapsar model a disk wind is thought to gen- 
erate the ^^Ni required to power the SN lightcurve 
(Macfadyen & Woosley 1999; Pruet et al. 2004). In 
the millisecond proto-magnetar model, the energetic 
wind shocks the material already processed by the 
supernova shock, perhaps generating the large in- 
ferred ^^Ni yields. Such a mechanism relies on tim- 
ing. We are currently investigating this scenario 
more fully. 

I thank Niccolo Bucciantini, Brian Metzger, 
Philip Chang, Eliot Quataert, and Jon Arons for 
stimulating conversations and collaboration. 
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